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1. Introduction 

Full counting statistics (FCS) is one of the most attractive and intellectu- 
ally involved concepts in quantum transport. It provides much information 
about charge transfer (all moments, including average current, noise, etc.) 
in a compact and elegant form. Albeit the study of full counting statistics 
has begun with some confusion. The first attempt to derive FCS [1] ex- 
ploited the text-book definition of quantum measurement: the probability 
of the outcome q of a measurement is given by (n\S(q — Q)\n), Q being an 
operator associated with the value measured and |n) the quantum state. 
However, the choice of Q made in [1] resulted in severe interpretation prob- 
lems. In [2] the authors revised their approach. The new method, that is 
commonly accepted now, invokes an extra degree of freedom, a detector. 
The quantum measurement paradigm is applied to the detector degree of 
freedom. 

There is a similarity between the method of [2] and two core approaches 
to quantum dynamics: the Keldysh technique [3] and the Feynman- Vernon 
formalism [4]. This similarity was not stressed in [2]. Later it has been 
noticed that the FCS can be evaluated with a straightforward modification 
of the Keldysh Green function technique [5] . This allows one to extend the 
studies of the FCS to many systems, some results being reviewed in other 
contributions to this book. This is good news. Seemingly bad news is that 
the method of [2] does not always give results that can be interpreted as 
probabilities of measurement outcomes. In the statistics of charge transfer, 
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the bad news becomes relevant for systems where gauge invariance is broken 
[6, 7]. If one tries to generalize FCS to arbitrary variables, these problems 
arise from the very beginning [8]. 

The present contribution consists of two parts. First we address the 
bad news. For this, we consider a rather general and abstract problem, the 
counting statistics of a general quantum mechanical variable being mea- 
sured by a linear detector. We find that in this case the FCS is physically 
meaningful and useful since it determines the quantum time evolution of 
the detector. This is despite the fact that it can not be interpreted as a 
probability distribution [8]. 

Thus encouraged, we investigate how far one can go with this approach. 
The abstract exercise with a linear detector proves to be very useful to 
describe a mesoscopic conductor embedded in a linear electric circuit. We 
reveal the relations between the FCS and the non-equilibrium Keldysh ac- 
tion that describes current and voltage fluctuations in electric circuits. We 
show how to evaluate the current and voltage counting statistics at any 
two contacts of an arbitrary circuit. This appears to be relevant for future 
experimental activities in the field of quantum noise and statistics [9]. 

The article is organized as follows. In section 2 we will pursue a route 
similar to that proposed in [1], defining a current statistics as the statistics 
of a charge operator. This route has not been explored before. It allows 
to understand the origin of the "negative probabilities" found in [6] purely 
in terms of conductor degrees of freedom. We follow an alternative route 
in section 3 and consider a system coupled to a static linear detector [8]. 
We show that the result of the detection depends in general on the initial 
state of the detector. To predict the result, one needs a function of two 
variables. We adopt then this function as the definition of the full counting 
statistics. We will identify situations where the two approaches give the 
same easy-to- interpret result. 

Then we turn to electric circuits. We explain in section 4 why a linear 
circuit can be viewed as a set of dynamical linear detectors and why the FCS 
expression is the building block of the Keldysh action that describes the 
quantum dynamics of the circuit. In section 5 we discuss the low-frequency 
limit of the action. In this limit, the action is local in time. This facilitates 
its evaluation. The general scheme is illustrated in section 6 with a simple 
model of a phase-coherent conductor in series with an external linear re- 
sistor. In this case, we are able to determine the full statistics of current 
and voltage fluctuations. These results were obtained in collaboration with 
C.W.J. Beenakker, that we gladly acknowledge. 
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2. Charge statistics without detector 

We start out by denning a statistics of time averages of electric current 
without specifying a measurement procedure, as it has been also done in 
[7]. Unlike in [7] we study directly the statistics of charge Q on one side 
of the cross-section of a conductor through which the current is measured. 
If Q is fixed at time 0, the statistics of Q(t) at time r corresponds to 
the statistics of the charge that has traversed the cross-section in the time 
interval [0,r]. Having the clear physical interpretation of the statistics of 
charge on one side of the cross-section, this definition will shed some light on 
interpretation problems encountered with other definitions, like "negative 
probabilities" and half-integer charge transfer [6, 7]. 

To specify Q we introduce, following [2], a smooth function / that di- 
vides the conductor into two parts: a left side, where /(x) < 0, and a 
right side, where /(x) > 0. We are interested in charge transfer through 
the boundary /(x) = between these two sides. The operator of electric 
charge to the right of the cross-section is Q = e9[f(x)] [where 8(x) = 
for x < and 0(x) = 1 for x > 0], e being the elementary charge. It is 
convenient to express the statistics of Q[r) via the corresponding moment 
generating function 



Xc(A) = E^T (QW) = ( eiAQ(T) ) = (e iHT e iXQ e- iHT ) • « 



Here H is the Hamiltonian of the current conductor and the average is 
taken over the initial density matrix of the conductor. 

One advantage of defining Xc as in Eq. (1) is that it is evidently associ- 
ated with a (positive) probability distribution, the distribution of measure- 
ment outcomes of an observable corresponding to the Hermitian operator 
Q(t). Besides, Xc predicts charge transfer in integer multiples of the elec- 
tron charge e for systems of well-localized non-interacting electrons, the 
result to be expected. Both very physical properties have been found to be 
violated by other definitions [6, 7]. One buys these advantages by defining 
a statistics that is only indirectly linked to the statistics of transmitted 
charge. That link requires the knowledge of the initial charge state of the 
conductor. 

Using that 



k 



we rewrite Eq. 



e iXQ/2 e A e -iXQ/2 = exp [ e ^Q/2^ e -UQ/2 ] 

(1) identically as 



(2) 



Xc (A) = (e lAQ / 2 exp [ ie - aQ ^He aQ / 2 r 




(3) 
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The charge operator Q commutes with all position operators x contained 
in the Hamiltonian H = iT(p,x). It does, however, not commute with 
the momentum operators p. As a consequence, momentum operators are 
transformed as 

e -iAQ/2 pe iAQ/2 = p _ | AV 0[^ X )] = pA ( 4) 

We define a new Hamiltonian H\ = H(p\,x) in the same way as it has 
been done in [2]. The generating function 

Xc (A) = (j*Q/* e iH-xT e -iH X T e i\Q/2\^ ^ 

takes then a form that is very similar to the x(^) found in [2] within the 
spin-^ current detection model . Tracing back the difference we write x(^) 
in terms of charge operators as 

X (A) = ^ e -<AO/2 e iAO(r) e -iAO/2^ (g) 

If the initial state of the conductor is an eigenstate of charge with eigen- 
value Qo, the two generating functions are identical up to the offset charge 
Qo, xW = e~* A< 3°Xc(A) . In this case both characteristic functions are 
associated with a positive probability 

P(q) = J d\e-* x X (\) (7) 

to transfer q charges during the measurement time. For generic systems 
of non-interacting electrons this probability has been found to be non- 
zero only at integer multiples of the electron charge, corresponding to the 
transfer of an integral number of electrons [2] . 

If the initial state is a superposition of eigenstates of charge Q, the 
two generating functions x an d Xc differ. For example, x ma Y seem to pre- 
dict the transfer of half the elementary charge for systems of well-localized 
non-interacting electrons, when Xc indicates integer charge transfer. This 
becomes evident when x is written in the form (6). It contains summands 
of the form Xmn = exp[— i\(m + n)/2] {m, a\ exp[iAQ(r)]|n, (3), \m,a) and 
\n, (3) being eigenstates of the charge with additional quantum numbers a, 
/?, Q\m, a) = m\m, a). While the matrix element in the expression for Xmn 
corresponds to integer charge transfer, the pre-factor suggests the trans- 
fer of half-integral charges when m + n is an odd number. This has been 
observed in [7]. 

For a Josephson junction at the cross-section /(x) = between two su- 
perconductors the situation is even more involved. One would like to choose 
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an initial state of well denned phase difference between the two sides of the 
junction in order to have a well defined current flowing. However, phase and 
charge are conjugated variables obeying Heisenberg's uncertainty principle 
and the dispersion of the initial charge Q in such a state is infinite. The 
charge distribution corresponding to Xc with an initial state of well-defined 
phase is therefore one of undetermined Q(t). It contains no information 
about the charge transfer. This is clearly undesirable. Because all the un- 
certainty of Q is already present in the initial state and has nothing to 
do with the transfer process, one might hope to be able to eliminate it in 
a meaningful way. Indeed, we have seen that this is accomplished for the 
initial charge offset in a state of well defined charge by the factors e~ tX ®/ 2 
in the definition of xW (6). One could hope that this works also for the 
initial charge spread in a phase eigenstate. Instead, x(ty then seems to 
imply "negative probabilities" [6]. 

A way to remedy this problem would be to calculate Xc in an initial state 
with a finite charge and phase dispersion. One could also try to find general 
relations between the charge statistics and the initial state of the conductor 
that would characterize the charge transfer process more generally. 

The alternative is to couple the (super-)conductor to a detector and 
interpret the detector read-off in terms of charge transfer. We will follow 
this route in the next section. It turns out that in a idealized detector model 
a function very similar to the generating function x eventually determines 
the final state of the detector provided the initial state is known. The 
charge transfer through the conductor is then characterized by its effect on 
the detector. 

3. Counting statistics with a static detector 

We turn now to the measurement process of time averaged quantities. In this 
section we focus on an idealized measuring device without its own dynamics, 
a static detector. Within this model we study the measurement of time 
averages dt A(t) of an arbitrary operator A. For electric currents, A = I, 
this will allow us to characterize the charge transfer and its statistics by 
a function similar to the characteristic function introduced in the previous 
section. The analysis of this function will allow us to trace the origin of the 
"negative probabilities" found in [6]. 

3.1. DETECTOR MODEL 

We employ a detector model that has first been used by John von Neu- 
mann in an analysis of the quantum measurement process. Following [10] 
e introduce a detector variable whose operator x commutes with all oper- 
ators of the system to be measured. Its Hamiltonian is q 2 /2m, where q is 
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the variable conjugated to x. The system measured shall be coupled to the 
position x of the detector during the time interval [0, r] and be decoupled 
adiabatically for earlier and later times. For this we introduce a smooth 
coupling function ot T {t) that takes the value 1 in the time interval [0, r] 
and is zero beyond the interval [ti,*2] (t\ < and £2 > t). The values for 
t± < t < and r < t < ti are chosen in a way that provides an adiabatic 
switching. The entire Hamiltonian reads then 

H{t) = H sys -a T {t)xA+^-. (8) 

The Heisenberg equation of motion for the detector momentum q 

q{t) = a T (t)A(t) (9) 

suggests, that the statistics of outcomes of measurements of the detector's 
momentum after having it uncoupled from the system corresponds to the 
statistics of the time average JJ" dt A(t) that we are interested in. 

The coupling term can be viewed as a disturbance of the system mea- 
sured by the detector. To minimize this disturbance, we would clearly like 
to concentrate the detector wave function around x = 0. The uncertainty 
principle forbids us, however, to localize it completely. Thereby we would 
loose all information about the detector momentum, which we intend to 
measure. This is a fundamental limitation imposed by quantum mechanics, 
and we are going to explore its consequences step by step. To discern it 
form a classical back action of the detector we take the limit of a static 
detector with m — > 00, such that x = and any classical back action is 
ruled out. 

3.2. APPROACH 

To predict the statistics of measurement outcomes in our detection model 
we need the reduced density matrix of the detector after the measurement, 
at t > t2- If there were no system to measure we could readily express 
it in the form of a path integral in the (double) variable x(t) over the 
exponential of the detector action. This is still possible in the presence of 
a system coupled to the detector [4]. The information about the system to 
be measured can be compressed into an extra factor in this path integral, 
the so-called influence functional. This makes the separation between the 
detector and the measured system explicit. To make contact with [2], we 
write the influence functional as an operator expression that involves system 
degrees of freedom only. We denote the initial detector density matrix (at 
t < t\) by p m (x + ,x~) and the final one (at t > £2, after having traced 
out the system's degrees of freedom) by pf(x + ,x~). R denotes the initial 
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density matrix of the system. The entire initial density matrix D is assumed 
to factorize, D = Rp m . 

We start out by inserting complete sets of states into the expression for 
the time development of the density matrix 



pf(x + ,x-) = Tr ( x +\l*e- i $ dt l. H '»- aT ® xA+ £\ D 

System 

y e'/t'i dt [ H ^~ a -^ xA+ 4;] \ x ~). (10) 



Here, T ( T ) denotes (inverse) time ordering. As the complete sets of states 
we choose product states of any complete set of states of the system and 
alternatingly complete sets of eigenstates of the position or the momentum 
operator of the detector. Those intermediate states allow us to replace the 
position and momentum operators in the time development exponentials 
by their eigenvalues. We can then do the integrals over the system states 
as well as the momentum integrals and arrive at the expression 

pf(x+,x-)= ! V[x + ] V[x~] p in [x + (ti),z-(ti)] e -iSvet([x+],[x-}) 

J X+(t2)=X+ X~(t2)=X- 

Tr Te~ i f*i dt [ Haya ~ aT ® x+ ® A ] R Te'^ dt[H ByB -a T (t)x-(t)A] ^-q 

System 



with the detector action 
We rewrite the expression as 



S Det ([x+], [*-]) = - jf* dt ^ [{x+f - (x-f] . (12) 



pf(x + ,x ) = J dxfdx 1 K(x + ,x ; xf , x 1 , r)p m (xf, x 1 ) (13) 
with the kernel 

K(x + ,x~;xf,Xi,T) = / V[x + ] P[aT] 

x+ (t2)=x+ ,x+ (t\)=x J [ X~ (t2)=X~ ,x~ (ti)=x± 

Z Sys ([a T x + ], [arx-]) e -«Det([x+],[x-]) (M) 
that contains the influence functional 

Z Sys ([ x + ], [x~]) = Tr Te^^i dt [ H ^-x + m] R f jjg dt [H syB - x - it)A] 

System 

(15) 
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Taking the limit of an infinite detector mass now, we find that <SDet i n Eq. 
(14) suppresses all fluctuations in the path integral. Therefore, the kernel 
K{x+ ■ X - X ^ - X . r) becomes local in position space, 

K(x + , x~ , xf, x± , t) = 5{x + — x~i) 5(x~ — x±) P(x + ,x~~,t) (16) 

with 

P(x + ,x~,t) = Tr re _i ^ 2d 4^ s y s_aT W x+A ]i?Te i ^ t i 2d< [ //sys_a!T ^ :E_j4 ]. 

System 

(17) 

It is constructive to rewrite now the density matrices in the Wigner 
representation 

P(x,q) = J ^e-^ z p(x+ Z -,x- Z -) (18) 

and define 

P(x,q,r) = j ^e-** P(x+ Z -,x- Z -,t). (19) 

With these definitions we find a compact relation that summarizes the 
results of this subsection: 

p f (x,q)= [ d qi P(x,q- qi ,T) p m (x, qi ). (20) 



3.3. INTERPRETATION OF THE FCS 

We adopt the relations (17), (19) and (20) as the definition of the FCS 
of the variable A. Let us see why. First let us suppose that we can treat 
the detector classically. Then the density matrix of the detector in the 
Wigner representation can be interpreted as a classical probability distri- 
bution Il(x, q) to be at a certain position x with momentum q. This allows 
for a classical interpretation of P(x, q, r) as the probability to have mea- 
sured q = Jq A(t). Indeed, one sees from (20) that pf(x,q) is obtained 
from p m by shifts in q, P(x, q, r) being the probability distribution of those 
shifts. 

In general, the density matrix in the Wigner representation cannot be 
interpreted as a probability to have a certain position and momentum since 
it is not positive. Concrete calculations given below illustrate that P(x, q, r) 
does not have to be positive either. Consequently, it cannot be interpreted 
as a probability distribution. Still it predicts the results of measurements 
according to Eq. (20). 



Full Counting Statistics in electric circuits 



9 



There is, however, an important case when the FCS can indeed be in- 
terpreted as a probability distribution. It is the case that P(x, q, r) does 
not depend on x, P(x,q,r) = P(q,r). Then, integrating Eq. (20) over x, 
we find 

Uf(q) = f dq' P(q-q',r)U m (q') (21) 

with H(q) = J dxp(x,q). Therefore, the FCS is in this special case the 
kernel that relates the probability distributions of the detector momentum 
before and after the measurement, II m (g) and n/ (q), to each other. Those 
distributions are positive and so is P(q,r). 

When studying the FCS of a stationary system and the measurement 
time r exceeds time scales associated with the system, the operator ex- 
pression in Eq. (17) can be seen as a product of terms corresponding to 
time intervals. Therefore in this limit of r — ► oo the dependence on the 
measuring time can be reconciled into 

P(x + ,x-,t) = e -£{x + ,*-y (22) 

where the expression in the exponent is supposed to be large. Then the 
integral (19) that defines the FCS can be done in the saddle point approx- 
imation. Defining the time average A = q/r, that is, A is the result of a 
measurement of JJ" A(t)dt/r, the FCS can be recast into the form 

P{x,A,t) =e- E(x ' A)T , (23) 

where £ is defined as the (complex) extremum with respect to (complex) 
z: 

£ = extr{£(x + ^-,x-^-)+iAz}. (24) 
z 2 2 

The average value of A and its variance (noise) can be expressed in terms 
of derivatives of £: 

d£(x + z/2,x - z/2) 



(A) = - lim 



z^o idz 
= ^(x-M/2x-z/2) 
z^o dz 2 

More generally, the quantity P(x + , x~ , r) is the generating function of mo- 
ments of q. It is interesting to note that in general this function may gener- 
ate a variety of moments that differ in the time order of operators involved, 
for instance, 

oM aN-M 

QZ = (- 1 )M^ ___ p (3; + x -, T) 



/' 

J 



dh...dt N {T{A(t 1 )...A(t M )}T{A{t M +i)-A{t N )}). (26) 
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The moments of (the not necessarily positive) P(0, q, r) are expressed through 
these moments and binomial coefficients, 

Q iN) - jdq q N P(0, q, r) = 2~ N £ (*) Q». (27) 

For A = I these moments correspond to the moments generated by the 
x(A) of the previous section, that is contained in P(x + , x~, r) as x(A) = 
P(A/2, — A/2, r) [compare Eqs. (5) and (6) in the semiclassical approxima- 
tion, where H\ = Hq — XI, to Eq. (17)]. Interpreting x(A) as the char- 
acteristic function corresponding to a probability distribution is therefore 
equivalent to the classical interpretation of P(x, q, r) discussed above. It is 
applicable to systems with an x- independent P(x,q,r). 



3.4. FCS OF A SYSTEM IN THE GROUND STATE 

To acquire a better understanding of the general relations obtained we 
consider now an important special case. We will assume that the system 
considered is in its ground state \g), so that its initial density matrix is R = 
\g)(g\- In this case the FCS is easily calculated. We have assumed that the 
coupling between the system and the detector is switched on adiabatically. 
Then the time development operators in (17) during the time interval t\ < 
t < adiabatically transfer the system from \g) into the ground state 
\g(x ± )} of the new Hamiltonian H sys — x^A . In the time interval < t < r 
the time evolution of the resulting state has then the form 

e -it(H sys -x±A) |^ x ±^ = g -itB(x±) | 5 ( x ±)). (28) 

Here, E{x ± ) are the energies corresponding to \g(x ± )). This gives the main 
contribution to the FCS if the measurement time is large and the phase 
acquired during the switching of the interaction can be neglected in com- 
parison with this contribution, 

P(x + ,X~,t) = e ~ir\E{x+)-E{x-)] m (2g) 

We now assume the function E(x) to be analytic and expand it in its Taylor 
series. We also re-scale q as above, A = q/r. We have then for the FCS 

P(x, A,r) = rjdz e~ izAr ■ e -ir[E'(x)z+E>"(x)z* /24+...] _ (30) 

First we observe that P(x, A, r) is a real function in this case, since the ex- 
ponent in (30) is anti-symmetric in z. A first requirement for being able to 
interpret it as a probability distribution is therefore fulfilled. However, the 
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same asymmetry assures that all even cumulative moments of A are iden- 
tically zero, whereas the odd ones need not. On one hand, since the second 
moment corresponds to the noise and the ground state cannot provide any, 
this makes sense. On the other hand, this would be impossible if P(0, A, r) 
were a positive probability distribution unless it had no dispersion at all. 

Belzig and Nazarov [6] encountered this situation analyzing the FCS 
of a super-conducting junction. In a certain limit the junction becomes a 
Josephson junction in its ground state. In this limit the interpretation of 
the FCS as a probability distribution does not work any longer. Fortunately 
enough, relation (20) allows us to interpret the results obtained. 

In the limit r — > oo of Eq. (30) all the terms involving higher derivatives 
of E{x) are negligible and we have 



According to the Hellman-Feynman theorem E'{x) = — {g(x)\A\g(x)). As 
one would expect, in this limit the measurement gives the expectation value 
of the operator A in a ground state of the system that is somewhat altered 
by its interaction with the detector at position x. Therefore the resulting 
dispersion of A will be determined by the initial quantum mechanical spread 
of the detector wave function. The error of the measurement stems from 
the interaction with the detector rather than from noise intrinsic in the 
measured system. 

3.5. FCS OF ELECTRIC CURRENT IN A NORMAL CONDUCTOR 

A complementary example is a normal conductor biased at finite voltage. 
This is a stationary non- equilibrium system far from being in its ground 
state. Here we do not go to microscopic details of the derivation. Our im- 
mediate aim is to make contact with the approaches of Refs. [2, 5]. We keep 
the original notations of the references wherever it is possible. 

The FCS of the current in a phase-coherent conductor is characterized 
by a set of transmission coefficients r ra [c. f. Eq. (37) of [2]]. At zero tem- 
perature and bias voltage V it reads 



lim P(x,A,t) = 5[A + E'(x)]. 



(31) 



P(x + ,X , T 



)=exp{^-\V\S[ie(. 



x 



.+ 



— X 




(32) 



with the function 




(33) 



n 



This expression depends on x + — x only, as a direct consequence of gauge 
invariance. Indeed, in every time development operator of Eq. (17) (with 



12 



MARKUS KINDERMANN AND YULI V. NAZAROV 



A = I) the coupling term is the coupling to a vector potential localized in 
a certain cross-section of the conductor. The gauge transform that shifts 
the phase of the wave functions by ex^ on the right side of that cross- 
section eliminates this coupling term. This transformation was explicitly 
implemented in [5]. Since there are two time development operators in the 
expression, the coupling terms cannot be eliminated simultaneously when 
x + 7^ x~. However, the gauge transform with the phase shift e(x + +x~)/2 
makes the coupling terms depending on x + — x~ only. 

Since P(x + , x~ , r) depends on x + — x~ only, the FCS P(x, q, r) is inde- 
pendent of x. As we have seen in section 3.3, under this condition the FCS 
can be interpreted as a probability distribution. As in section 2 we con- 
clude that for the statistics of the current of non-interacting electrons the 
characteristic function %(A) = P(X/2, — A/2, r) is associated with a positive 
probability distribution. 

Superconductivity breaks gauge invariance, thus making such an inter- 
pretation impossible. 

4. Electric Circuits as General Linear Detectors 

The model used in the previous section may seem rather abstract and 
unrealistic. To make contact to the "real world", we notice that 

— the time derivative of the detector momentum q is related to the cur- 
rent through the mesoscopic system 

— the velocity x enters the Hamiltonian in the same way as the voltage 
applied to the conductor. 

Next we adopt the following definition of the " real world" : the only quanti- 
ties measured are electric voltages and currents between nodes of an electric 
circuit. The most adequate description of the quantum mechanics of the 
system would thus contain these variables only. This description is hardly 
possible to achieve within a Hamiltonian formalism, since the latter can 
contain neither dissipation nor retardation. It is the Feynman- Vernon for- 
malism [4] that allows to formulate quantum mechanics in the form of an 
action that contains only necessary variables. This action may be derived 
from the Hamiltonian formalism by tracing out irrelevant degrees of free- 
dom. But this is precisely what we have done in the previous sections! 
The conclusion is that the above results can be used to arrive at an ade- 
quate formulation of the quantum dynamics of electric circuits that contain 
mesoscopic conductors. This formulation clarifies the notion of a detector. 

We formulate the problem as follows. Suppose we know the FCS of 
a mesoscopic conductor. When measured, it is embedded in an electric 
circuit. Generally speaking, one does not measure the voltage or current 
directly at the mesoscopic conductor but rather somewhere else in the cir- 
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a) b) 




Figure 1. Electric circuit as a general linear detector. The black box symbolizes the 
electronic system to be measured. It is embedded in a linear environment (shaded box). 
Either voltage (a) or current (b) is measured between the output contacts 1 and 1' of 
the circuit. 



cuit. We cast this in Fig. 1 in circuit theory terms. Either a voltage or a 
current measurement is performed at the output contacts 1 and 1'. The 
mesoscopic conductor (black box) is connected to the input contacts 2 and 
2'. The shaded box presents the electromagnetic environment of the meso- 
scopic conductor. It is supposed to be a piece of a linear circuit and thus 
can be characterized separately from the "black box" by three (frequency- 
dependent) response functions. The question is what is the FCS of such a 
measurement. 

We answer it by extending our simple detection model of the previous 
section to a set of dynamical detectors. We model the "environment", in- 
cluding the detectors, by a set of additional degrees of freedom. This way 
we study classical back action effects of the detectors. These effects will 
be more important than the subtle "quantum" influence of the detector 
on the measured system discussed in the previous section. The detectors, 
representing the environment, shall obey linear equations of motion. This 
restriction to linear circuits is compatible with most experimental situa- 
tions. Detector non-linearities are avoided in many experiments because of 
undesired effects like the mixing of different frequencies. 

4.1. THE MODEL OF A LINEAR ELECTROMAGNETIC ENVIRONMENT 

Linear systems can be modeled by a set of non- interacting bosons [11]. The 
Hamiltonian reads 

^linear = ^ ^m a Jn a m ~ ^ ( 34 ) 
m i 

where the operators Hi are the physical quantities of interest and conju- 
gated to generalized forces hi. These operators are linear combinations of 

Hi = J2 c m a ™ + h - c ~ ( 35 ) 
m 
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A proper choice of Q m , chl models arbitrary response functions and enables 
one to formulate the quantum dynamics of the circuit in terms of quantities 
of interest only [11]. We are dealing with two pairs of conjugated variables, 
(Xi,Ii) an d (x2>-^2)) where Xj(t) = J drVj{r). One variable of each pair 
can be chosen as generalized force and the remaining one as operator. Al- 
though different choices lead to different Hamiltonians, the concrete choice 
is just a matter of convenience. Our primary choice is to treat I 2 and xi 
as generalized forces. The response functions relate I\ and V 2 to those two 
variables: 

Ji(w) = y(w)Fi(w) + tf(w)j 2 M 

V 2 (u) = Z(u)I 2 (u) - K^Viiu). (36) 

More symmetric choices express either currents in terms of voltages via 
admittances Yi,Y 2 , Y12, 

h(uj) = Y 1 (u)V 1 (u)+Y 12 (o J )V 2 {u;) 

I 2 (uj) = Y 2 (u)V 2 {u) +Y 12 {u)V 1 (u), (37) 

or voltages in terms of currents via impedances Z\,Z 2 , Z\ 2 , 

ViH = Z^hiu) + Z 12 {u)I 2 {u) 

V 2 {uj) = Z 2 (u)I 2 (u) + Z 12 (u)h(u). (38) 

There are obvious relations between these response coefficients: Z = 1/Y 2 ,Y = 
(YiY-2 - Y&)/Y 2 , K = Y 12 /Y 2 , (Z u Z 2 , Z l2 ) = (Y 2 , Y u -Y 12 )/(Y 1 Y 2 - Y&). 

Here we assume a passive circuit in thermal equilibrium. The response 
functions then satisfy Onsager symmetry relations and the fluctuation- 
dissipation theorem relates the response functions to fluctuations of the 
corresponding variables. 

4.2. GENERAL RELATION 

We derive now the general relation that determines the full counting statis- 
tics of the output I\ (t) of the linear detector described above provided it is 
coupled to the "black box" to be measured. We will follow the lines of sec- 
tion 3.2. However, now we assume the detector to be in a state of thermal 
equilibrium at temperature T. In addition, instead of evaluating the final 
density matrix of the detector p* , we analyze the partition functional Zj 
that generates moments of the read-off environment variable I\ . We define 
this functional as 
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It generates moments of outcomes of measurements of I\ (t) at any instant 
of time. 

It is advantageous to write the functional in its dependence on the 
combinations Vf = {d/dt){xf + Xi~)/2 and x\ = {xt — Xi)/^i that we 
collect into a "vector" xi = (^i d >Xi)- The "classical" field Vf accounts for 
a (time-dependent) voltage source between the contacts 1 and 1'. Moments 
of I\{t) are generated by differentiation of Zi with respect to the anti- 
symmetric, sometimes called "quantum", field: 

([AteM-Ml) = z^.--zs4i^KJL^ (40) 

In general, these moments depend on Vf 1 . 

Hi is the Hamiltonian of the circuit for a current measurement. It reads 

Hi = linear + #Sys - X2I ■ (41) 

-^linear governs the bosonic detector degrees of freedom and i^Sys acts on 
the electronic degrees of freedom of the "black box" . The third term couples 
the electric current to the detector degree of freedom X 2{t) = J dt' V 2 {t'). 
The latter is thus the analogue of x in section 3. 

Both X2 and I\ are linear combinations of boson creation/annihilation 
operators, 



X2 



E c (x) a i r (x)* t j Ji) , J/)*J (42) 



We rewrite now Zj as a path integral in detector variables, like we have 
done to derive Eq. (11). The integration variables now are a,m\t), ± cor- 
responding to the two parts of the Keldysh contour. Operators X2,h arc 
replaced by 

X? = £ Wo® + h-c. l{ ±] = £ 4?c4 ±) + h.c, (43) 

m m 

the sign depending on the part of the contour they reside on. To proceed, we 
introduce two extra vector variables into the path integral: X2 = (V^^xt) 
and q\ = (^f ><?2)' by inserting the identity 

1 ~ h[V?]V[xl\X[5{2V?-xi-X2)mxl-xi + X2) 
J t 

~ J V[V 2 d ]V[x q 2 }V[I?}V[q«] 

exp |i J dt [If {2x1 ~xt + X2) ~ q q 2 m d -xt- X2 )]} • (44) 
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This allows us to split the integrand into two factors. One factor is a trace 
over "black box" variables. It depends on X2(t) only and is in fact the in- 
fluence functional i?Sys introduced in section 3.2, Eq. (15) (with A = I). 
It generalizes our definition of the FCS of the "black box" biased by the 
voltage V d (t) to time dependent arguments. The other factor is a quadratic 
form in the a m (t). Its Gaussian character enables us to perform the inte- 
grations over the a m , yielding 

Zj[xi] = J P[x 2 M£>K l{5cnv(Kl] ^^ (45) 
<5env = \% (xiYXi + 2XiKq 2 + (46) 
Scoup = 2 J dt [ql{t)V?{t) - If{t)x q 2 {t)] • (47) 

The environmental part of the action is expressed in terms of response 
functions of our primary choice. The 2x2 matrices Y, K in the time domain 
are integral kernels depending on the time difference only. In frequency 
representation they read 

Y(u,\ ( ° \ (AR\ 

K > \Y(u) -2za;coth(^)Rey(a;) ) ' { ' 

K( \ ( -K*{u) \ ,.qs 

K[UJ) ~ \K(u) 2cjcoth(^)ImK(a;) )' { Uj 

The matrix Z is of the same form as Y. We can further simplify this relation 
by integrating over qi . The result acquires the more symmetric form 

Zi[xi] = J V[Ue- lS - (mm) Z^ s [^ (50) 
^ [xi^iXi + 2x1^12X2 + £2^2X2] , (51) 



where the matrices Y\,Y<2,Yi2 are defined as in Eq. (48). This is the de- 
sired general relation: it expresses the FCS of electric currents through a 
mesoscopic conductor that is embedded in a linear electric circuit. It is im- 
portant that all information about the mesoscopic conductor enters through 
its FCS at voltage bias, i%ys- 

Let us discuss the relation in some detail. First, let us replace the general 
four-pole electric circuit by a single two-pole resistor Z(uj) (see Fig. 2). We 
do this by choosing Y\ = Y 2 = —Yi 2 = Y = l/Z. The resistor and the 
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O 




Z(co) 



V 2 X 



Figure 2. A simple example of a linear environment: a resistor Z(u>) in series with the 
"black box". The circuit is biased with a voltage source V. The Quantum dynamics of 
the system is formulated in terms of fluctuating fields V"2,X2- 



"black box" enter the expression then on equal footing, 

ZliXl] = J V [X2\ ^Resistor [X2 ~ Xll^Sysfe], (52) 
^Resistor K] = exp | -1 j ^X^xj , (53) 

where Y is defined as in Eq. (48). -Zrteistor defines the Gaussian FCS of the 
linear voltage-biased two-pole resistor. Eq. (52) expresses a simple concate- 
nation rule for the FCS of compound circuits: the total FCS is a convolution 
of the FCS's of the individual elements. The resulting FCS is presented as 
an integral over the fluctuating phase X2 defined in one node of the cir- 
cuit. This is in agreement with known results about Keldysh actions of 
electric circuits with Josephson junctions [12]. Moreover, the rule can be 
easily generalized to more complicated circuits, for example, circuits with 
two mesoscopic conductors in series. 

It is interesting to note that the integration over the fluctuating phase in 
fact implements a constraint, namely current conservation, on the quantum 
motion of the compound system. The current through the "black box" 
equals the input current of the detector: 

<5(/Sys " h) - Il^SysW " itWHWsysV) - I{ (t)) ~ 
t t 

J I%exp j* J dt[xt(t)(I S ys(t)-lt(t))-X 2 (t)(ISys(t)-I{(t))]Y 
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the last equality holding for operators. Since the two currents are repre- 
sented by operators of a different nature (fermionic for Is ys and bosonic 
for Ii), the constraint would be difficult to handle within a Hamiltonian 
formulation. This illustrates the usefulness of the Keldysh action approach 
in this context. 

So far we have disregarded the constant factors that arise on various 
stages of the derivation from the variable changes in the path integrals 
and/or subsequent Gaussian integrations. This is not a matter of our care- 
lessness but rather indicates a problem with the normalization in the ap- 
proach in use. The normalization of the final relation (50) is known from the 
fact that Zj(x) = 1 at x q = P er definition of the generating functional. 
In the path integral representation the correct normalization is assured by 
the causality structure of the action and integrals over the whole range 
of frequencies contribute [13]. One can not reproduce the correct factors 
within the low-frequency approximations and the saddle-point treatment 
of the path integrals employed in this article. The problem is not related 
to the embedding of the mesoscopic conductor and persists also to com- 
pound linear circuits. One may draw analogies with known, general and 
long-standing problems related to the quantum mechanics of constrained 
motion [14]. To our knowledge it is always possible to correct for wrong 
factors that arise in calculations by "normalizing" the results to Zj(x) = 1 
at x 9 = 0- F° r an general relations given in this article we assume that 
the correct factor is included into the definition of the metrics in function 
space, that is, into the definition of V[x2\- 

4.3. VOLTAGE MEASUREMENT AND PSEUDO-PROBABILITIES 

The evaluation of the FCS of electric current through a mesoscopic conduc- 
tor with voltage bias is a rather straightforward task. One proceeds within 
the Hamiltonian formulation in terms of electronic operators subjected to 
an external voltage. The conjugated problem is the FCS of voltage fluc- 
tuations under the condition of current bias. The absence of an obvious 
Hamiltonian formulation for current bias makes this problem less straight- 
forward. 

Albeit the problem can be solved in a general way within the approach 
of the previous subsection. Let us consider the generating functional of 
voltage fluctuations between 1 and 1' (Fig. 1 b) defined by 

Z v \q{\ = (je^ dt [Hv+^(t)Vi]^ e -ifdt {Hy+q+iW}^ _ (54) 

Similarly to Eq. (40), derivatives of this functional with respect to q q (t) 
give the moments of voltage fluctuations. We now repeat the derivation of 
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the previous subsection for this functional. The answer reads: 
Zy\q\\ = J V[q 2 ]V[x 2 ]e~ 1 ^^^ 

Serw = J ^ [qiZiqi + 2q{Z 12 q 2 + fa^fc] , 



■'coup 



= 2 I dt[ql{t)V 2 d {t)-I${t)xl{t) 



(55) 
(56) 
(57) 



The matrices Z\ 2 \ 2 {lo) are defined analogously to (48) with the impedances 
(38). 

We now choose the linear detector in such a way that the voltage fluc- 
tuations between 1 and 1' are the same as between 2 and 2' and the con- 
ditions of current bias are met. Elementary circuit theory tells us that this 
is achieved in the limit Z\ 2 —> —Z\ = —Z 2 — > —00. In the expression (55) 
this corresponds to setting q\ = q 2 and <S cnv = 0- 

This brings us to a remarkable conclusion: for any conductor, the gen- 
erating functionals for voltage noise at current bias I cl (t) and for current 
noise at voltage bias V cl (t) are related by a functional Fourier transform, 



(58) 



Therefore, the functionals Zy and are in fact just different forms of the 
same object. There is some uncertainty in this definition stemming from 
the problem mentioned in the previous section. Zy has to be "normalized" 
such that Zy = 1 at q q (t) = 0. We assume the normalization factor to be 
included into the definition of V. 

This simple relation between Zy and Zi suggests to define the functional 



H[v],[i}) 



I 



v[qQ] e 2ifdt q *V Zy 

V[ X q ] e 2i f dtx91 Z! 



I 

qi 
V 

X" 



(59) 



that depends on "classical" variables only. We show with the help of Eq. 
(40) that P has the following properties: 



and 



(V(t 1 )...V(t m )) 



<I(*i).../(t m )> 



_ JV[V}V(t l )...V(t m )P([V],[I}) 
lit)' JV[V]P([V],[I]) 

jv[i)i(t l )...i(t m )P{[vui}) 

v(t)~ JV[I]P([V],[I]) 



(60) 



(61) 
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Eqs. (60) and (61) resemble the properties of a probability density and yet 
P is not a probability. First, it does not have to be positive, just as the 
P(x,q,r) discussed in section 3 needs not. Second, P is dimensionless in 
contrast to a true probability density of either voltage or current fluctua- 
tions. For these reasons we refer to P as a 'pseudo-probability'. 

P([V], [/]) is a characteristics of a two-pole conductor and takes a simple 
form for a linear circuit. In this case P is positive and depends only on the 
impedance of the conductor Z(u) and its temperature T, 

From Eq. (60) one derives the voltage correlator 

(V(tf)= J ^Re^)coth^, (63) 

that conforms to the fluctuation-dissipation theorem. 

The general relation (50) can be rewritten in terms of pseudo-probabilities. 
In this case, it will express the pseudo-probability of current/voltage fluc- 
tuations between 1 and 1' i\([V], [I]) in terms of the Ps ys of the "black 
box" and a four-pole pseudo-probability Pyi that depends on two currents 
and voltages, characterizing the linear part of the circuit: 



Pl([Vl], lh}) = J P[^2p[/2]iWl], [h\; [V 2 ], [h])Psy S ([V2], [h])- (64) 

We cast now also the result Eq. (52) for the "black box" in series with 
a linear resistor Z(uS) in terms of pseudo-probabilities. The answer is ex- 
pressed in terms of the pseudo-probability of the resistor defined by Eq. 
(62): 

A an n = / v[v 2 \p z ([v - v 2 \, wp^m, m 

The above relations are transparent if one interprets them in terms of 
classical probabilities. They show that the probability of a certain cur- 
rent/voltage fluctuation is composed of a probability of fluctuations in the 
"black box" and of a probability of fluctuations in the linear detector, 
voltage and current satisfying circuit theory rules. Yet the relations are 
quantum mechanical ones and are written for pseudo-probabilities. They 
contain all necessary information about quantum properties of the system 
under consideration. 

To summarize, fluctuations in a mesoscopic system that is embedded 
in a linear environment can be characterized by three statements. First, 
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voltages and currents are related by circuit theory rules. For instance, for 
the circuit of Fig. 2 the current fluctuation SI in the "black box" results in 
the current fluctuation 51/(1 + ZG) in the whole circuit, G being the linear 
conductance of the " black box" . This factor is rather trivial from a classical 
point of view. However, in the context of quantum mechanical detection, it 
constitutes the main source of detector back-action. This "classical" back- 
action could be disregarded for the static detector considered in section 3. 
Second, the linear part of the circuit produces its own Gaussian quantum 
and thermal fluctuations. Also these fluctuations contribute to the FCS of 
the compound system. Third, the fluctuations in the linear environment 
affect the fluctuations produced by the " black box" . This effect is the least 
trivial and we will explore its consequences in section 6. 

5. Low-frequency limit 

We have shown in the previous section that the FCS of a mesoscopic conduc- 
tor in a general linear environment can be expressed by the relation (50). 
This relation involves path integration and defines a quantum field the- 
ory with the field x*2(t). This field theory is non- linear, the non-linearities 
coming from the mesoscopic conductor. The latter is sometimes underes- 
timated, since the I — V curve of the conductor may be perfectly linear. 
However, even in this case the noise and higher order correlators in general 
do depend on the voltage, which provides the non-linearity. 

Such non-linearities make the general problem already hardly tractable. 
To complicate the situation further, the most important part of the action, 
ln^sys, is only known for stationary %. In this case it is given by the rela- 
tion (32). Microscopic Keldysh Green function techniques can be used to 
evaluate this functional at non-stationary x- The answer has, however, a 
complicated functional dependence on the argument \. 

Fortunately enough, there is a simple low-frequency limit where the 
quantum field theory becomes tractable. This limit applies to many exper- 
imental situations. In this limit we consider only quasi-stationary realiza- 
tions of the fields. The action for these realizations is local in time, which 
allows for easy path integration. 

Before turning to the concrete formulation, let us discuss the time scales 
that occur in our problem and that may determine the non-locality of the 
action. There may be a time scale, the i?C-time, characterizing the fre- 
quency dependence of the detector response functions Yi,2,i2- This scale 
depends on the system layout and can be easily changed. From the exper- 
imental point of view it is convenient to choose the response functions to 
be constant in a wide frequency interval. This allows us to disregard this 
time scale. Another time scale, tq, has quantum mechanical origin. The 
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environment action (51) becomes non-local at frequencies uj ~ /cbT, where 
the crossover from classical to quantum noise occurs. For the mesoscopic 
conductor a similar scale may arise from the energy eV related to single 
electron transfer through the conductor. Therefore tq ~ 7i/max{eV, fc^T}. 
A third time scale is the typical time interval tj between electron transfers 
through the mesoscopic conductor, e/I. Comparing the latter two scales, 
we find that tj ~ (e 2 //iG)rg, G being the conductance of the "black box". 
Generally, we expect the low-frequency limit to be valid at time scales ex- 
ceeding both 17 and tq. These time scales are typically short in comparison 
with the time resolution of the measurement electronics. Consequently most 
electric measurements are performed in the low- frequency limit. 
In this limit, the general relation (50) takes the following form: 

Zi([Vi],[xi]) = Jv[V 2 ]V[ X 2]exp{- J dt[£ cnv (V l (t),xi(t);V 2 (t), X 2(t)) 

+S Sys (V 2 (t), X 2(t))}}; (66) 

where 

WVi , Xi ; V 2 , X2) = ~H(ViYi + V 2 Y l2 ) X i + (V 2 Y 2 + V 1 Y 12 ) X 2] 

+ AT(Y lX l + Y 2 x 2 2 + 2Y 12 xiX2), 
£ Sys = -\n Z Sys (V,x)/r. 

In the last equation, Zs ys is evaluated for stationary V, x- Thus its logarithm 
is proportional to the measurement time r (see Eq. (32)). For the resistor 
in series (Fig. 2) £ env reduces to: 

£c„v = £z(Vi -V 2 , X i- X2); £z(V, X ) = + 4T X ) X /Z. (67) 

As expected in the low-frequency limit, the noise of the detector is just 
Johnson-Nyquist white noise. 

A simple estimation shows that the action is not only local in time. It 
is also big for the time intervals in question. This allows us to implement 
a semiclassical approximation, that is, to evaluate the path integrals in 
saddle-point approximation. Indeed, the typical values of X , manifesting 
the discreteness of charge, are of the order of 1/e. Typical values of the 
action are then estimated as S ~ Ixt — I/er ~ (r/r/). By definition of 
the low-frequency limit, t/ti S> 1. 

We now come back to the central problem addressed in this article: the 
distribution of currents averaged over the time interval r. Taking the path 
integrals in saddle-point approximation, we arrive at a simple expression 
for the FCS of the "black box" embedded in an environment circuit: 

Zi{V uX i) = exp[-£ 7 (Vi,xi)T], 

£i{Vi,Xi) = exti{£ cnv (V l ,xi;V 2 , X 2) + £sys(V 2 ,X2)}. (68) 

V2,X2 
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The corresponding FCS of voltage fluctuations at given current I can 
be obtained directly from Eq. (68) by using Eq. (58), 

Zv(I,q) = exp[-£V(I,g)r], 

S v (I,q) = exti{£ I (V, X )-2i(I X -qV)}. (69) 

V,x 

When integrating Eq. (58) we take care of the normalization to assure 
that Zy(I,0) = 1. A similar integration determines the pseudo-probability: 

P(I,V) = exp(-f(V,/)r), 

S(I,V) = extr{£ I (V, X )-2iI X }. (70) 

Finally, we give here expressions for real probabilities to measure the voltage 
V at a given current I, Py(I), and to measure the current 7 at a given 
voltage V, Pi(V). By virtue of the relations (60) and (61) we obtain 



Pv(I) \ _ [Z 
Pi(V) / " V 2t 




P(I,V). (71) 



Remarkably, these probabilities are in fact the same with exponential ac- 
curacy. They differ in normalization factors only, 



Py(I) _ dP PS 

Pr(F) V d 2 £dv 2 ' { ' 

The relations (69-72) are valid for the statistics of the outputs of the com- 
pound circuit as well as for the "black box" biased by either a voltage or a 
current source. 

The FCS in the low-frequency limit can thus be formulated in terms of 
stationary functions £i,£v,£(I-,V). These functions resemble thermody- 
namic potentials. As those potentials, they are related by Legendre trans- 
forms. However, the relations presented concern non-equilibrium systems. 
They provide an effective minimum principle for fluctuations in such sys- 
tems. 



5.1. WHERE ARE THE CHARGING EFFECTS? 

The electron transport through mesoscopic circuits is known to be affected 
by charging or Coulomb blockade effects. They are due to electron-electron 
interactions at the mesoscopic scale. The simplest case of a tunnel junction 
embedded into an electromagnetic environment is described in [15] in some 
detail. Although these charging effects are not the topic of this article, we 
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find it important to note that they are present in our formalism. Indeed, 
the Keldysh action in use is a generalization of the action introduced in 
[12] to describe Coulomb blockade effects in Josephson junctions. 

How do charging effects manifest themselves in the low-frequency limit 
considered? 

The way to understand this is to invoke the field-theoretical concept of 
renormalization. We have obtained the low-frequency action just by sub- 
stituting slow realizations of the fluctuating field xi into the original ac- 
tion. The renormalization procedure gives a more accurate way to obtain 
the low-frequency action. In this procedure, the field is separated into a 
slow and a fast part. The path integration over the fast variables is per- 
formed. Since fast and slow variables are coupled, this alters the functional 
of the slow variables, it is "renormalized" . For the problem under consider- 
ation, this renormalization introduces charging effects. Since the environ- 
ment action is quadratic, it is not renormalized. The renormalization can 
be thus ascribed to £sys- The typical time scale of relevant fast modes is 
tq. The importance of the charging effects is governed by the value of the 
environment impedance at the corresponding frequencies, Z(1/tq)[15]. If 
Z{\/tq) <C fi'/e 2 , the renormalization is small and charging effects can be 
disregarded. In the opposite case, charging effects are important and the 
renormalized action £ differs significantly from the bare one. 

We thus conclude that charging effects do not change the low-frequency 
relations given in the previous subsection. Their only effect is a renormal- 
ization of £ in comparison to its bare value at voltage bias. This renormal- 
ization depends on the environment response functions at the frequency 
scale 1/tq. 

6. Mesoscopic conductor in a macroscopic circuit 

In this section, based on Ref. [9], we illustrate the general relations obtained 
above. We do this by considering a specific example: a phase-coherent con- 
ductor in series with an impedance Z (Fig. 2) [9]. Further we focus on 
the shot noise limit eV S> ksT, and disregard the Johnson- Nyquist noise 
produced by the linear resistor. First we address the FCS of current fluctu- 
ations in the low-frequency limit discussed above. We will thus make use of 
the relations given in section 5. The assumption Z(lo ~ eV) <^iti/e 2 allows 
us to disregard charging effects without any restrictions on Z(0). 

To employ our scheme, we first need the concrete expression for the 
FCS of the voltage-biased conductor, i?sys- In the stationary case, this has 
been derived in [2]. The mesoscopic conductor is completely characterized 
by its transmission eigenvalues T n [16] that are assumed not to depend on 
energy. In our notations, we have (\eV\ S> fc^T) 
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2 S ys(V, X ) = exp{-£ Sy sT}, £ Sys = - ^ I V\ S[2ie X sign(F)] (73) 

with the dimensionless S of Eq. (33). 

Since we disregard the noise of the resistor, we have the environment 
action £ Z (V, \) = -2W Z~ x x- By virtue of Eq. (68), the FCS of the circuit 
is defined by 

£i(V, X ) = extr \-2iZ-\V 2 - V)( X 2 - x) ~ ^S[2ie X 2 sign(V 2 )] ) . 

(74) 

Minimization with respect to V 2 fixes X2, 

2iZ' l (x-X2) = ^S(2ie X 2), (75) 

and makes further minimization with respect to this variable redundant. 
We immediately arrive at 

£i(V,x) = -^S[2ie X 2(x)} (76) 

where X2(x) is implicitly defined by Eq. (75). In the limit of a vanishing 
external impedance Z — > 0, we have X2 — > X so that £/(V, x) — * ^Sys, as it 
should. 

The coefficients of a series expansion of £j in x give the cumulants 
of the transmitted charge ((Q p )), whereas the coefficients of £s JS are the 
same cumulants {{Q p ))o in the limit Z — > 0. Comparing the two series, we 
obtain relations between these cumulants. The linear terms imply (Q) = 
(l+ZG) -1 (Q)o, G being the conductance of the mesoscopic conductor. This 
is a consequence of Ohm's law for the average current I that is reduced by 
a factor 1 + ZG by the resistor Z in series. The Langevin approach of Ref. 
[17] predicts the same re-scaling for the fluctuations of the current. Indeed, 
in second order we find that {{Q 2 )) = (1 + Z 2 G)" 3 ((Q 2 )) , in agreement with 
this prediction. 

However, for higher order cumulants we find terms that are not con- 
sistent with this re-scaling hypothesis. For instance, the third cumulant 

llnH fglo 3ZG (((Q 2 )) ) 2 

W " (1 + ZG) 4 (l + ZGf (Q) ' [ ' 

Although the first term on the right-hand-side does have the scaling form 
conjectured, the second term does not. This is generic for p > 3: {(Q p )) = 
(1 + ZG)~ p ~ 1 ((Q p ))q plus a non-linear (multinomial) function of lower cu- 
mulants. This mixing in of lower order cumulants is a consequence of the 
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non-linear feedback mechanism discussed in the end of section 4.3. To give 
another example, the forth cumulant reads 

«Q 4 » - {m ° 



(1 + Z 2 Gf 



10ZG ((Q 2 ))o((Q 3 ))o 15Z 2 G 2 ((Q 2 ))g 



(l + Z 2 Gf (Q) (l + Z 2 Gy (Q)l 

We now turn to a one-channel conductor with transmission probability T. 
In this case, the distribution of the integer transmitted charge q = Ire for 
voltage bias V is known to take the simple binomial form [2] (4> = eVr/2-K) 



V.) - C) 



T q (l - T)*- q . (79) 



The distribution dual to this is that of the flux <p = (e/h) J dtV(t) that 
is accumulated at the conductor under current bias /. We make use of 
to the relation (70) that determines these probabilities. For a mesoscopic 
conductor in the shot noise limit, the corresponding 8(1, V) reads 

S(I, V) = extr | ^S(2ie X ) - 2i/ X } (80) 
Expanding this near its point of extremum, we find that 

di 2 d 2 s i 2 



d 2 £ dV 2 V 2 



(81) 



for any functional form of S. Employing Eq. (72) we readily establish that 
the distribution function of flux P q ((f)) is given by (here we set q = Ir/e) 



Pq (<P) = ^(q)=^_ l ^(l-T)t>-\ 



(82) 



which is known as the Pascal distribution. 

We complement this derivation by a simple interpretation of the Pascal 
distribution (82) of voltage fluctuations. The binomial distribution (79) of 
transferred charge can be interpreted [2] in gambling terms: it gives the 
probability to win (transfer an electron) q times in 4>q game slots given a 
winning chance of T. In the voltage biased conductor, electrons leave the 
leads with a frequency eV/h and 4>o = (e/h) JJ" dtV(t) is the number of 
electrons ("game slots") that try to pass the barrier in the detection time 
window. 
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The current bias changes the rules of the game. The gambling now ends 
when precisely qo attempts have been successful. This requires a fluctuat- 
ing number <p of game slots (electrons that hit the barrier). Indeed, the 
Pascal distribution quantifies the probability distribution of the number of 
independent trials that one needs to achieve a given number of successes in 
a Bernoulli experiment. 
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